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ABSTRACT 

Although  organized  transportation  systems  are  very  old  indeeds, 

active  research  into  their  operations  is  relatively  new0   Some  of  the 
problems  encountered  in  transportation  research  are  discussed  along 
with  methods  for  their  solution,.  The  Hitchcock  distribution  problem 
is  discussed  and  a  sample  problem  solved  using  three  different  methods. 
The  usefulness  of  graph  theory  in  studying  transportation  problems  is 
pointed  out.   Examples  are  then  given  of  its  use  in  a  maximal  network 
flow  problem  and  in  the  problem  of  minimizing  equipment  requirements 
to  meet  a  fixed  schedule.,   Finallys  some  miscellaneous  techniques  are 
discussed. 

The  writer  wishes  to  express  his  appreciation  for  the  guidance  and 
encouragement  given  him  by  Professor  Thomas  E»  Oberbeck  of  the  U„  S„ 
Naval  Postgraduate  Sehools  for  the  helpful  suggestions  of  the  second 
reader^  Professor  Ho  C„  Ayerss  and  for  the  careful  clerical  assistance 
provided  by  Mrs„  Norma  A„  Stevens „ 
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INTRODUCTION 

In  the  past  decade  research  Into  transport  At  lor,  problems  has   been 
greatly  intensified,,   Much  of  this  xesearch  has  involved  the  application 
of  the  principles  and  techniques  from  a  wide  variety  of  mathematical 
disciplines,,   In  this  thesis  we  shall  summarize  some  of  the  existing 
techniques  and  give  examples  of  their  applicability,, 

Most  of  the  research  we  have  found  reported  has  as  its  object  the 
optimization  of  some  feature  of  a  transportation  activity 0   Some  examples 
are;  minimizing  costSj,  relieving  congestions,  minimizing  equipment  re- 
quirement ss  establishing  economical  maintenance  schedules;,  minimizing 
loading  and  unloading  times „   Usually  the  methods  reported  are  directed 
at  solving  a  problem  of  a  particular  transportation  activity  such  as  a 
railroad  or  an  airline,,   But  much  of  the  research  has  been  sufficiently 
general  to  be  widely  applicable,,   The  tendency  to  generalize  models  is 
prevalent .   For 'example 9  a  researcher  may  seek  to  minimize  the  number 
of  "units"  required  to  maintain  a  fixed  schedule.   "'Units1"'  may  tnen  be 
interpreted  by  a  railroad  as  boxcars  or  trains;  by  an  airline 8  as  air- 
planes; and  by  a  shipping  company;,  as  tankers  or  freighters  „ 

Transportation  systems  of  any  appreciable  size  are  extremely 
complex  entities  and  may  defy  analysis  in  many  respects.   Researchers 
turn  to  mathematical  models  of  the  systems  they  wish  to  analyze  and 
build  into  the  model  the  degree  of  detail  or  generality  necessary  to 
give  significant  results.  This  practice  has  resulted  in  the  application 
of  a  great  many  mathematical  techniques  to  problems  in  transportation. 
The  techniques  and  disciplines  which  have  been  used  include ;  linear  pro- 
gramming,, transhipment  and  assignment  ti     ues8  non-linear  and  dynamic 
programming,,  queuing a  or  waiting  line  the©rys  graph  theory^  and  the  Monte 
Carlo  method.  1 


In  our  research  for  this  paper  we  ha  ye  found  a  great  deal  of  mater- 
ial on  transportation  type  problems  in  the  following  publications;  The 
Jou'rna  1  of  the  Operations  Research  Society  of  America,,  Management  Science,, 
Naval  Research  Logistics  Quarterly^  and  others  listed  in  bibliographies 
at  the  end  of  each  chapter.   Groups  we  have  found  to  be  particularly 
active  in  transportation  research  include:'   The  Rand  Corporation;  lanta 
Monica j,  California ;  the  George  Washington  University  Logistics  Research 
Project;  The  Cargo-Handling  Research  Project  at  the  University  of  Calif- 
ornia at  Los  Angeles;  The  Management  Sciences  Research  Group s  Purdue 
University;  to  name  only  a  few. 

Professor  F„  Harary  of  the  University  of  Michigan  has  done  a  consider- 
able amount  of  research  in  graph  theory9  which  is  a  directly  supporting 
discipline.   He  is  co-author  of  a  new  book  on  this  subject  which  is  due 
to  appear  in  1960,   This  book  will  be  the-  first  full  treatment  of  this 
subject  in  English, 

As  an  indication  of  the  direction  in  which  current  research  is  tend- 
ing,,, we  point  to  Rand  report  R-351  which  givess  in  abstract  form,,  the  pro- 
ceedings of  "The  Rand  Symposium  on  Mathematical  Programming".   In  his 
address  to  the  opening  session  of  this  symposium^  G,  B,  Dantzig  cited  the 
five  following  theoretical  areas  as  those  in  which  future  developments 
are  expected ; 


Special  Structures 

(2)  Discrete"  Programme  r.« 

(3)  Network  Theory 

(A)  Non- linear  Programming 
(5) 


It  is  anticipated  that  transportation  theory  will  be  one  of  a  great 
many  beneficiaries  of  developments  in  these  areas. 


In  arranging  various  types  of  transports^  ?s 

for  the  ensuing  discussion  we  ware,  of  necessity,  a  bit  arbitrary.  We 

have  presented  in  detail  only  a  few  of  the  more  widely  used  techniques 

and  have  made  reference  to  others „  We  have  tried  to  present  a  reason- 
ably complete  bibliography  at  the  end  of  each  chapter,, 

In  Chapter  I  we  have  stated  and  discussed  some  methods  for  the 
solution  of  the  Hitchcock  distribution  problem,,  both  with  and  without 
capacity  limitations,,   These  are  essentially  linear  programming  pro- 
blems.  In  Chapter  II  we  have  indicated  the  applicability  of  linear  graph 
theory  to  transportation  problems  and  general  network  problems*   In 
Chapters  III  and  IV  we  have  we  have  discussed  in  some  detail  two  problems 
that  illustrate  the  applicability  of  graph  theory  to  transportation 
problems.   In  Chapter  V  we  have  cited  several  examples  of  other  techniques 
and  the  problems  to  which  they  have  been  applied.   Specifically^  we  have 
discussed  queuing  theorys  simulation  and  straightforward  engineering 
analysis  including  the  difficult  and  time  consuming  task  of  data  gather- 
ing. 


CHAPTER  1 
THE  TRANSPORTATION  PROBLEK  -   I   EAR  PROGRAMMING 

1„   The  Hitchcock  Distribution  Problem,, 

The  problem  to  be  discussed  in  this  chapter  was  originally  stated 

in  1941  by  F„  L»  Hitchcock  Qj  and  was  farther  discussed  and  elaborated 

by  To  Co  Koopmans  [_2  J  „   The  problem  may  be  seated  as  follows.  Certain 

amounts,,   a   (i  =  Ij,,,^),  of  a  homogeneous  product  are  available  at 

each  of  m  origins^  and  certain  amounts  b,   (j  =  i,,oliliu)  are  required 

J 

at  each  of  n  destinations.^   The  costj,   c.  .9    of  shipping  a  unit  amount 
from  the  i   origin  to  the  j   destination  is  known  for  all  i  and  j„ 
We  will  further  stipulate   ^'.a,  -   «Z.b£«   The  problem  is  to  select  the 

i      j 

amounts  x.   to  be  shipped  form  each  origin  to  each  destination  so  that 

the  total  shipping  cost  is  &   minimum »   The  amounts  x„   must  be  non- 

ij 

negative „      xs  ■£■  0„ 
ij 

Algebraically  j,   we   state  the  problem  in  ich*;    following   forms 
Given?      (1)      Two   sets  of  constants     a  i  ■■  lSooa,m 

and  b  „        j  -  18 „ . .  sis 


2"ai  "^bj» 

(2)  A  matrix  of  constants  C  -   f  c.  .  J  i  «  1-, 


,m 


Find,;   A  set  of 


x»  .  ^  02  satisfy: 


^.  x.  .     -     b  J  ;r- 

c       ij             j  "' 

Js   i»- 

J>    x  „  „      -     a  „ 
-S-     ij              i 

i-    l,oo.,ni 

3 

ers  in  brackets  refer  to  references  <  !  bed  in  the  Bibliogra 


^c)     ^—  ^      c  .  „   K.   -  minimum 

»    »     i  t!   ij 

Linear  programming  problems  which  may  be  stated  in  the  foregoing 
form  have  come  to  be  called  "the  transportation  problem"'  even  though 
they  do  not  all  deal  with  transportation,,   Examples  of  such  problems  are 
the  personnel  assignment  problem^  the  contract  -  award  problem,,   (_3j 
f/4]  the  Traveling- Salesman  problem  [loj  8  and  the  Caterer  problem  Q.1J  , 

Go  Bo  Dantzig  [5 J  was  first  to  formulate  this  problem  as  a  linear  - 
programming  problem^  and  he  gave  an  algorithm  for  its  solution  based  on 
his  Simplex  method  [6 J  „   The  transportation  algorithm  is  considerably 
easier  to  apply  than  the  general  simplex  method  In  that  it  does  not  re- 
quire the  inversion  of  matrices „   The  algorithm  is  discussed  In  more 
detail  and  a  sample  problem  worked  out  In  Appendix  Ac 
2„   The  Assignment  Technique » 

As  was  noted  above 9  a  special  case  of  this  "transportation  pro- 
blem" is  the  n  x  n  assignment  problem,,   In  this  case  all  a  -  b„  -  18 
and  all  x,  ,  are  either  one  or  zero,,   The  problem  mav  be  stated  as 
follows;   determine  the  optimum  assignment  of  n  persons  to  n  jobss 
given  a  matrix  of  constants  [r„  „J  which  specify,,  in  appropriate  unitss 
the  worth  or  performance  rating  of  the  i    individual  In  the  j   job„ 
The  most  straightforward  approach  is  to  consider  each  of  the  n,  possible 
arrangements „   Buts  for  only  moderately  large  n„  this  becomes  a  hope- 
less tasko  What  is  required  is  a  procedure  for  finding  an  optimum  assign- 
ment (there  may  be  more  than  one)  with  a  reasonable  amount  of  effort « 
Such  a  procedure  was  provided  by  H„  Wo  Kuhn  |_3J  in  his  "Hungarian  Method", 
■°s  algorithm  has  been  modified  and  extended  to  the   m  x  n 


transportation  problem  by  L„  R*  Ford,  Jr„  and  Dc  R,  Fulkerson  f?J  , 

and  by  Jo  Munkres  [_8J  „   The  algorithms  so  developed  may  be  described 

as  procedures  that  make  maximal  use  of  minimal  cost  routes.   Each  assigns 

quotas  only  to  the  least  expensive  feasible  routes  until  the  assignment 

is  complete o   The  resulting  assignment  is  one  which  gives  minimum  costo 

It  is  not  necessarily  unique.   These  algorithms  are  also,  discussed  more 

fully  in  Appendix  A  and  a  problem  is  solved  using  each  in  turn, 

3,   The  Transhipment  Technique, 

An  interesting  variation  on  the  m  x  n  transportation  problem  (  m 

sources,,  n  destinations)  is  the  transhipment  problem  described  by  Alex 

Orden  [93   In  the  original  transportation  problem  each  point  acted  as 

a  shipper  only  or  as  a  receiver  onlys  and  the  routes  were  considered  to 

be  direct  from  each  shipper  to  each  receiver „   In  his  extension  of  the 

problems  Orden  has  permitted  shipments  to  go  via  any  sequence  of  points 

rather  than  being  restricted  to  direct  connections  from  one  of  the  origins 

to  one  of  the  destinations.   His  method  of  solution  involves  converting 

the  extended  problem  to  the  original  problem  by  &   rather  simple  device. 

Each  point  is  treated  as  a  pair  of  points 8  one  acting  as  a  shipper  and 

one  as  a  receiver,   The  unit  cost  of  shipment  from  a  point  considered 

as  a  shipper  to  the  same  point  considered  as  a  receiver  is  set  equal  to 

zero,   Ail  other  c  „   (i9  j  ~  190„0,,  m  +  n«   i£?j'»  are  assumed  known 

ij 

and  greater  than  zero.   For  computational  purposes  large  stockpiles  are 
assumed  to  exist  at  each  shipping  point.  The  excess  of  stockpiles  over 
actual  amounts  shipped  appear  as  shipments  from  a  point  to  itself  at 
zero  cost  and  are  removed  from  the  final  result.   With  this  mechanism^, 
the  m  x  n  problem  with  transhipment  is  converted  to  an  M  x  M  problem 
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at  trans  >ck  problem) s   where  M  -  m  +  o„  The 

problem  i  oni  ands  in  the  non-    raerate  casei(   2M 

-  1  routes  have  positive  quotas  assigned.   It       a  chat  in  minimum  - 

cost  solutions  all  x. ,  are  positi     bull  th       removed  leaving  M  ■ 


1  =  m+n-1  actual 

4,   Integer  Solutions, 

It  is  worth  'noting  here  that  if  the  Hitchcock  problem  is  posed  in 
integers  (i.e.  each  a   and  b„   is  a      Lve  Integer)  then  there  is 
at  least  one  minimizing  solution  In  which  each  x.   is  either  a  positive 

J*  J 

integer  or  zero  l_4j  ,   In  this  case  the  whole  ptoblem  can  be  solved  in 
integers  by  introducing  s      feasible  solutions  consisting  of  integers 
and  then  changing  these  only  by  integral  amounts.   Very  recently,,  Gomory 

j_13j  has  developed  a  "method  of  integer  forms'"   for  solving  a  general 
linear  programming  problem  with  the  additional  constraint  that  the 
variables  in  the  solution  be  integers,      tzig  [14J  has  included  this 
procedure  In  the  term  "discrete  programming"  and  forsees  in  its  further 
development  promise  :\:    solving  "all  kind?  of  combinatorial ,  non- linear 9 
non-convex  problems, ,, areas  "where  classical  mathematics  has  been  weakest"  = 
5„   The  Capacitated  Problem, 

In  the  classic  Hitchcock  problem  the  upper  bound  on  the  amount  of 
commodity  assigned  to  any  route  is  determined  by  the  aaiounts  available  at 
the  sources^   a  'Sj   and  the  amounts  required  at  the  destinations,,   b/s. 
The  routes  themselves  are  not  regarded  as  limited  in  their  capacities. 
One  can  imagine  practical  situations  wherein  the  actual  capacities  of  the 
routes  may  be  limiting  factors.   For  examples,  suppose  It  is  desired  to 

*See  Appendix  A  for  a  discussion  icy. 


transport  certain  known  »m  several  sources  to 
several  destinations  -.-      ■yi    or  .    available  certain 
railway  routes  with  specified  capacities;  I  •  i          -  of  the  com- 
modity shipped  on  route  i-js  in  the  specif  Line  intervals  cannot  ex- 
ceed the  capacity  of  that  route „ 

Our  introduction  c               I  the  problem,, 
at  least  to  a  limited  degree,,   The  natural  specification  of  capacity  is 
a  rate*   Consider  a  prot  iem   sim  cock  problem^  except  sup- 
pose the  a  " s  represent  quanr  ipped  in  a  certain  time  in- 
1 

terval9  say  a  month,,  and  the  b  „°  s  are  amounts  required  at  the  destina- 

J 

tions  per  month„   Further8  let  us  spe        capacity  of  each  route  i-j 

in  units  of  commodity  per  month,,  i  ths  so  -  called  capacitated 

Hitchcock  problem,,   This  problem  is  simiiai  itated  problem,, 

except  that  we  envision  it  as  bet  ■    applicable  within  time 

•  intervals  oi        . A   lengthy  Basicall  -,  it  Ls  static  in  that  we  have 

introduced  no  truly  dynamic  varial         as  instantaneous  flow  rates 

or  transit  times.  And  we  have  suppc  coi  ;tants-a„  °  s8  b„8ss 

1  J 

and  the  capacities-   do  not   change 
We  state   the  problem  formal: 


[a)      Two   sets  of  cons    i 
li 


as  i  2   Ijooojin 


»  D   o   o  o  5)  lltl 


that       >    a,      =      >_   b .   =  I 

(b)  A  matrix     C     =  fc.  .  J 

L    ij  J 

(c)  anothi  rix     L     =     d 


where  c    represents  the  capacity  of  ..id 

d  ,   represent  8  the 

amount  from  source  i  to  destination  j. 
Find  X  ■   f"  x .  .  1 

such  that 

(1)   5 
C2)   5" 


3    iJ     1 

(3)  0  £  xtj  £  ctj 

(4)  ^.  ^L  x.  *d,  ,   =  minimum 

Whereas  the  incapacitated  problem  always  has  a  feasible  solution  |_4J  : 
the  capacitated  problem  may  not,   The  problem  is  feasible  if,  and  only  if, 
the  maximal  flow  of  the  associated  network  is  at  least  equal  to  the  total 
amount  to  be  shipped,,   K„   A  discussion  of  networks  and  their  maximal 
flows  will,  follow  in  Chapters  II  and  III„ 

Ford  and  Fulkerson  have  also  provided  a  method  for  solving  the  capa- 
citated problem  in  their  paper,   "A  Primal  -  Dual  Algorithm  for  the  Capa- 
citated Hitchcock  Problem""  QsJ  ,    As  their  title  implies.,  these  authors 
define  a  dual  problem  involving  dual  variables  q/      £   ,,   and  "j^ 
These  are  used  to  distinguish  certain  routes,  and,,  in  terms  of  them,  a 
restricted  primal  problem  is  solved.   From  this  solution,  new  dual  vari- 
ables are  defined  and  the  restricted  primal  problem  is  again  solved,   this 
time  using  the  new  dual  variables.   In  each  solution  of  the  restricted 
primal  problem.,  a  set  of  positive  x„  .   is  determined,,   When  the  sum  of 


these  x, .   is  equal  to  the  total  amount  to  be  shipped^   K,  the  x; .   so 
determined  provide  a  solution  to  the  original  problem. 

During  each  iteration;  in  determining  the  solution  to  the  restricted 


9 


primal  problem,  the  authors  call  upon  a  special  flow  algorithm  developed 
by  them  in  a  previous  paper  (j.6J  »      The  mechan       I  «e  solution  are 
similar  to  those  described  for  the  incapacitated  problem  and  will  not 
be  described  in  further  detail  here. 

In  a  subsequent  paper  (_17J  s  Fulkerson  showed  that  an  upper  bound 
on  the  number  of  iterations  or  "labelxngs""  required  to  solve  the  trans- 
portation problem  could  be  determined.,  These  upper  bounds  were?   (a) 
for  an  n  x  n  optimal  assignment  problem,,   n{n+l);   (b)   for  the  tran- 
shipment problem,   R(n+1);   (c)   for  a  Hitchcock  -  Koopmans  problem 
K  Cmin(m,n)+l3  . 

The  literature  contains  numerous  other  examples  of  applications  of 
linear  and  non- linear  programming  techniques  to  transportation  problems. 
Several  such  papers  are  cited  in  the  bibliography  without  special  reference 
to  them  in  the  texto  These  are  intended  to  be  representative^  rather  than 
exhaustive,  to  attempt  an  exhaustive  listing  would  be  presumptuous s  indeed. 
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CHAPTER  11 
TRANSPORTATION  NETWORKS  AND  GRAPH  THEORY 

1,   Discussion. 

In  the  stud/  of  transportation  systems  it  bie  to 

utilize  network  diagrams.   The  usual  approach  is  to  let  the  nodes  of  the 
network  represent  origins^  destinations,  or  junction  points,,  and  to  con- 
nect these  nodes  with  lines  representing  the  available  routes.  The  des- 
cription, representation^  and  use  of  such  diagrams  has  been  greatly  facili- 
tated by  the  use  of  the  mathematical  theory  of  linear  graphs.   An  early 
example  of  the  use  of  graph  theory  in  c     t ion  with  transportation 
problems  was  in  a  paper  by  T.  C  Koopmans  and  S.  Reitei  [  ]  j       :ied 
in  1951 „   Much  of  the  work  had  been  done  earlier  during  World  War  II 
by  the  former  author. 

The  time  -  lag  in  the  application  of  graph  theory  to  transportation 

problems  is  emphasized  by  quoting  from  Che  above  mentioned  paper 

"The  cultrual  lag  of  economic  thought  in  the  applicati     i  .mathe- 
matical methods  is  strikingly  illustrated  by  the  fact  that  linear 
graphs  are  making  their  entrance  into  transportation  theory  just 
about  a  century  after  they  were  first  studied  in  relation  to  elec- 
trical networks,  although  organized  transportation  systems  are 
much  older  than  the  study  w£  electricity." 

Perhaps  some  of  this  delay  is  due  to  the  fact  that8  u?    /ery  recent- 
ly s,  graph  theory  has  been  somewhat  neglected  by  mathematicians  themselves., 
Few  articles  on  the  subject  appear  in  mathematical  journals. 

Until  1958y  there  was  only  one  book  devoted  entirely  t    raph  theory 
-  Do  Konig"s  Theorie  der  Endlichen  und  Unend lichen  Graphen  {^2 J  3    first 
published  in  Leipzig  in  1936  and  reprinted  by  Chelsea  in  New  York  in  1950. 
The  second  book  on  the  subject s  by  Claude  Serge  ^3  J   B  was  published  in 
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France  In  1958. 

Professor  F,  Harary,  of  the  Umi  is  writ  tea 

quite  extensively  on  graph  theory  ands  with  Rc  Z,   Norman,  is  writing  a 
book  due  to  appear  in  1960c   [5 J 

Because  of  its  wide  applicability  and  comparative  simplicity,  the 
two-terminal     network  is  used  in  many  transportation  problems.   Normal- 
ly, flow  in  such  a  network  is  unidirectional  from  a  source  to  a  sink. 

This  theory  has  been  especially  useful  In  determining  the  maximal 
flow  which  a  network  will  sustain.   Examples  of  its  use  in  this  con- 
nection are  included  in  Chapter  III, 

The  similarities  between  communications  and  transportation  net- 
work problems  have  been  recognized  and  pointed  out  by  R,  E.  Kalaba  and 
M,  L,  Juncosa  £4  J  ,  among  others.   Ibis  Similarity  Is  not  limited  to 
their  common  tie  with  graph  theory  but  extends  also  to  queuing  theory ? 
simulation  devices,  linear  programming,  dynamic  programming  and  Boolean 
algebra, 
2,   Representations, 

We  shall  exclude  from  our  discussion  graphs  containing  arcs  which 
join  individual  nodes  to  themselves, 

A  natural  representation  of  a  network  is  a  diagram  in  which  points 
represent  nodes.,  and  lines  represent  the  routes  between  them.   Elsewhere, 
(Appendix  C)  we  have  defined  a  network  t©  be  a  (connected)  graph  G,  to- 
gether with  the  capacities  of  its  individual  arcs.   We  shall  illustrate 
here  a  one-origin,  one  terminal  planar  network  with  source  a  and  sink  b„ 

1,   For  definitions  of  this  and  other  technical  terms,  see  Appendix  C 
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The  graph  is  not  oriented, 


a 


F: 

This  network       also  1  a  symmetric  matrix 

L  =\JL       \     in   which  (     is  the  capaci  c  ij  if  the  arc  ij 

ij  ij 

is  included  in  the  graph s  and  is  zero  otherwise 0   'for  our  example  the 

applicable  matrix  would  be  the         ■  Except  on  the  main  diagonal. 

os  are  not  written  but  are  implied  by  a  blank  space. 
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In  Chapters  III  and  IV  we  shall  use  Che  linkage  or         i  of 
graphs.   For  represents",  rams  such 

as  figure  2-1  rather  tl     he  matrix  repre  2-2s  will 

be  used. 

It  is  desirable  to  emphasize  that  we  have  pointed  out  only  one  field 
for  the  application  of  graph  theory.   In  ac  I  Lt   i  Co  their  usefullness  In 
the  study  of  transportation,,  communications 9  an         al  networks s 
graphs  have  been  used  In  many  other  fields.   They  are  useful  in  describ- 
ing Markov  chains,,  for  example.   As  another  example.  Pa      r  Harary 
has  made  use  of  them  in  the  study  of  the  behawv  a  1  -sciences  £o3  9  C?J  , 
Much  of  his  work  in  graph  theory  has  been  don  L ion  with  the 

Research  Center  for  Group  Dynamics,  University  of  Mi  higan  and  8 
Social  Science  Research,,  Bell  Telephone  Laboratory 
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CHAPTER  III 
MAXIMAL  NETWORK.  FLOW 

1,   General 

As  we  mentioned  in  Chapter  IIS  transportation  systems  are  often  re- 
presented by  network  models.,   The  nodes  may  represent  simple  transfer 
points  such  as  the  pick  up  point  for  the  hook  in  a  cargo-handling  opera- 
tion; orj,  they  may  represent  highly  complex  entities  such  as  a  seaport 
or  a  complete  railway  operating  division,,   What  the  nodes  represent  de- 
pends upon  the  model  under  consideration  and  the  degree  of  generality 
or  detail  under  which  it  is  being  studied.   These  nodes  are  joined  by 
arcs  which  may  have  an  associated  traffic-carrying  capacity.   This 
collection  of  points  and  arcs  with  their  capacities  will  be  referred 
to  as  a  network  -  (see  Appendix  C) , 

One  question  of  interest  is  that  of  determining  the  maximal 
steady- state  flow  through  such  a  capacitated  network.  We  shall  des- 
cribe in  this  chapter  two  procedures  for  determining  maximal  flow.,  one 
a  flooding  technique  due  to  Alexander  W,  Boldyreff  [ij  an<^  t^®  other 
the  minimal  cut  theorem  due  to  L.  R„  Ford,  Jr,  8  and  D,  R,  Fulkerson  |_2j  , 
2C   The  Flooding  Tjachn_i_que 

Boldyreff  applied  his  technique  directly  to  traffic  through  a  rail- 
road network^,  and  was  careful  to  claim  no  generality  for  it.   He  character- 
ized it  as  empirical  in  nature,,  and  left  as  an  open  question  the  possibil- 
ity of  wider  use  of  the  flooding  technique.   It  appears  feasible  to  apply 
it  in  any  situation  where  a  network  model  will  fairly  describe  the  actual 
situation. 

The  method  employed  can  be  applied  to  networks  with  any  number  of 


origins  emd   terminals;  howevers  Coi  clarity  the  one-origin,  one-ter- 
minal case  Is  described. 

Certain  subsets  of  points  ate  distir.  gins,  the 

terminals j,  and  junction  points  (see  Appendix  C),   We  will  discuss  flow 

in  a  railway  network  assuming  'capacities,  in  trains  per  day„   By  "steady- 
state"  we  mean  that  the  flow  rate  through  each  arcs  in  trains  per  day>, 
is  constant  from  day  to  day„   The  steady- state  is  further  clarified  in 
these  initial  assumptions; 

lo  At  the  origins  only  loaded-  trains  leave  and  only  empty  t tains 

arrive, 
2a   At  the  terminals  only  loaded  trains  arrive,  and  only  empty 

-ains  leave „ 
3,  At  each  junction  point  the  number  of  loaded  (empty)  trains 
arriving  is  equal  to  the  number  of  loaded  (empty)  trains 
leaving  each  day„ 
40   The  number  of  loaded  trains  Leaving  all  the  origins  is  equal 

to  the  number  arriving  at  all  the        Is  each  day„ 
50   The  number  of  empty  trains  leaving  all  the  terminals  is 
equal  to  the  number  arriving  at  all  th<    Lgins  each  day. 

The  last  two  statements  are  deducible  £r  i   -  first  three  and  are  stat- 
ed explicitly  only  for  emphasis,, 

The  problem  is  to  find  a  method  of  assigning  the  steady- state 
flow  of  traffic  to  each  arc  of  the  network.,  not  exceeding  the  capacity 
of  the  arcs  which  will  maximize  the  total  flow  of  loaded  (empty)  trains 
from  the  origins  (terminals)  to  the  terminals  (ori  s  :isfy 


the  five  assumpt     stated  *b: 

In  the  problem  as  e       the  st«  iists  of  simultaneous 

movement  through  the  network  of  loaded  I  origins  to  the 

terminals  and  empty  trains  in  the  opposite         i,   These  two  move- 
ments must  necessarily  be  at  the  same  rate.   By  imagining  a  network 
with  arc  capacities  just  one-half  of  those  of  the  original  network,  we 
can  look  only  at  the  flow  of  loaded  trains  frc     gins  to  terminals  and 
reduce  the  problem  to  the  unidirectional  case.   By  unidirectional  flow  we 
mean  that  the  flow  in  any  one  arc  is  in  one  direction  only.   In  this  dis- 
cussion the  flow  considered  will  be  generally  directed  from  the  origins 
to  the  terminals „ 

For  this  simplified  problem  the  steady-state  conditions  are  con- 
tained in  the  following  statements; 

1.   At  the  origins  all  trains  leav  Lve. 

20   At  the  terminals,,  all  trains  an.     none  depax 

3o   At  each  junction  point  the  number  of  trains  coming  in  is 
equal  to  the  number  going  out. 

4„   The  number  of  trains  leaving  t  equal  the  number 

arriving  at  all  termina". 

5,   The  traffic  flow  through  each  arc  cannot  exceed  the  capacity 
of  that  arc 
By  defining  two  networks  to  be  equivalent  if     naximal  flows  through 
them  are  equal,-,  we  can  readily  reduce  a  network  with  many  origins  and 
many  terminals  to  an  equivalent  network  with  a  single  origin  and  a 
single  terminal. 

'We  shall  describe  the  flooding  techniq    utilizing  a  single-origin, 

2D 


single- terminal  network,,   The  procedure  is  as  follows; 

Starting  at  the  origin,,  assign  sufficient  traffic  flows  to  all  arcs 
leaving  the  origin  to  saturate  them.   This  gives  the  maximum  number  of 
trains  arriving  at  junction  points  one  arc  removed  from  the  the  origin., 

View  this  set  of  junction  points  as  new  origins  andy  starting  with 
the  one  subject  to  the  greatest  capacity  constraint s    schedule  trains 
whenever  possible  in  the  following  order; 

lo    "Forward"  -  to  new  junction  points  through  the  outgoing  arcs 

2,  "Laterally"  -  to  other  points  of  the  set, 

3.  "Bottlenecked0  -  if  trains  are  left  over  after  steps  1  and  2; 
i,e,  if  all  outgoing  and  lateral  arcs  are  saturated* 

At  such  a  set  of  junction  points  there  may  be  arbitrary  decisions 
regarding  which  trains  to  "bottleneck".   The  guiding  principle  xs   to 
move  forward  as  many  trains  as  possible8  and  to  maintain  the  greatest 
flexibility  for  the  remaining  network. 

Continue  the  above  procedure  until  the  scheduled  flow  covers  the 
complete  network  and  reaches  the  terminal „ 

Eliminate  bottlenecks  by  returning  all  excess  trains  to  the  origin. 
The  validity  of  the  solution  can  be  checked  by  inspection.   If  a  maximal 
flow  has  been  achieved9  there  will  be  no  continuous  unsaturated  path 
extending  from  the  origin  to  the  terminal,,   (That  is  there  will  be  no 
chain  connecting  origin  and  terminal  that  does  not  contain  at  least  one 
saturated  arc).   If  this  criterion  is  not  met8  the  inspection  will  have 
revealed  the  unsaturated  chains.   Flow  in  these  chains  can  be  increased, 
giving  maximal  flow.   As  will  be  seen,  the  observation  concerning  satur- 
ated arcs  is  a  point  of  similarity  between  this  procedure  and  the 
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ii8<     ia 

-    »aper8  Boldyi  lot    simplifying 

complicated  networks,   H     >,  as  he  points  owt9  the  networks  to  which 
these  procedures  can  be  applied  are  rather  snore  the  exception  than  the 
rule,-,  and  the  applications  are  tedious,   *'"-  e  author  recommends  a  straight- 
forward application  of  the  flooding  techniq  i 

A  sample  problem  is  solved  fey  this  procedure  in  Appendix  B,   The  same 
problem  is  also  solved  using  Ford  and  Full    n°8  minimal  cut  theorem, 
3,  The  Minimal  Cut  Theorem, 

This  method  was  applied  to  networks8  in  general,,  not  necessarily 
rail  networks.  Although  we  will  not  give  hers  a  proof  -  minimal 
cut  theorem  8  we  will  need  some  de£:.  .  a   order  to  discuss  its  ap- 

plications.  See  Appendix  C, 

A  gjrapjh,,  Gs  is  a  finite^  one-dim-     :  al  complex^  composed  of 

vertices  a8b8c8o,»8e  and  arcs   <X.  (ab)s  yS  (ac)s...s  cT(ce), 
An  are  CK    (ab)  joins  its  end  vertices  aj,b;  it  passes  through  no  other 
vertices  of  68  and  intersects  other  arcs  onl^  in  vertices, 
A  chain  is  a  set  of  distinct  arcs  i       which  can  be  arranged  as  ck   (ab)? 

&   (bc)8,,08  j  0  (gh  where  the  vertices  s  fc  c8,,,8h  are  distinct,  A 
chain  does  not  intersect  itself;  i  t ices  a  and  h.   Each 

arc  in  G  has  associated  with  it  a  j  bex  called  its  capacity. 

The  graph  G8  together  with  th    ipa  -.dividual       is 

called  a  network.   Two  vertices   i     ar^  distinguished  as  the  source  and 
b  the  sink,   A  chain  flow  from  a  to  to  i       iple   (Csk)  composed  of  a 


1.   For  a  proof  of  the  thoreim9  see  the     ad  article  listed  in  the 
bibliography  at  t    i  •  -i!.  of  this  cl  ■  ■  r. 
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chain  joining  a  and  b  and  a  non-negativ     ber  k  repress. u     the  flow 
along  C   form  sour      sink,  a  flow  in  a  network  la  a  collection  of 
chain  flows  which  has  the  property  thaft  the  sum  of  the  nunil  -    f  all 

chain  flows  that,  contain  any  arc  is  no  greater  than  the  capacity  of  that 
arc.   If  equality  ho  Ids  s  the  arc  is  said  to  be  saturated  by  the  flow,, 
A  chain  is  saturated  with  respect  to  a  flew  if  it  contains  a  saturated 
arc   The  value  of  a  flow  is  this  sum  of  all  the  chain  flows  which  compose 
it.  A  disconnecting  set  Ds  is  a  set  ©£  arcs  which  has  the  property  that 
every  chain  joining  a  and  b  meets  the  collection.   A  disconnecting  set, 
no  proper  subset  of  which  is  disconnecting.,  is  called  a  cut . 

The  value  of  a  disconnecting  set  v(D»  Is  the  sum  of  the  capacities 

of  its  individual  members.   Thus  ad-  g  set  of  minimal 

value  is  automatically  a  cut. 

We  state  the  minimum  cut  theorem. 

Theorem  I,      The  maximal  flow  value  otota        n  a  network  is  the 
minimum  of  v(D)   taken  ever  all  d£  s  Do 

This  theorem  is  so  intuitively  appealing  that  it  hardly    ams  to 

require  proof  at  all.   However^  as  examj  -    n  of  r«      :<s  [,2 J  will 

show,,  the  proof         involved,,  arri  for  certain  in  c  example 

see  (f.3j  )  it  does  not  apply.   The  inapplicability  to  networks  with 

several  sources      >rresponding  sinks  wa^  Ford  and 

Fulkerson  [2j  ■.  # 

e  actual  computation  procedure  for  determining  tl      .i   of  the 

minimal  cut  is  based  on  a  corollary  to  tl  i  minima  L         em  and  an 

al  theorem  which  will  be  stated  wit 

Corollary?   Let  A  be  a  colli  ■  res  of  a  network  N 
wh:             cut  of  N  in  If  N"  is  a  net- 

work obtained  from  N  by  ad  £   «*ach 

a re  o f  As  then  cap    ° 
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eorem  2.   If  the  graph  of  Gs  together  with  arc  abs  Is  a 
planar  grapn  CO  ,  there  existfi  a  chain  joining  a  and  b  which 
meets  each  cist  of  N  precisely  once0 

Let  T  be  the  chain  joining  a  and  b  wl       I  upmost  in  \h      T  has  the 

property  of  theorem  2.   Impose  as  large  a  chain  flow  as  possible  (T^k)  on 

this  chain,,  thereby  saturating  one  of  its  arcs. 

By  the  corollary,  subtract  k  from  the  capacity  of  each  arc  of  To 
Remove  the  previously  saturated  arc  whose  capacity  is  now  zero.   Continue 
this  procedure*   Eventually  the  graph  di  -    seta  and  a  maximal  flow  has 
been  established. 

For  an  Illustrative  example 9  see  Appendix  B. 

In  problems  of  sea- transport,  the  most  restrictive  capacity  limitations 
are  quite  often  the  cargo- hand ling  capacities  of  ports  at  which  ships  must 
load  and/or  unload^  rather  than  on  sea  routes  themselve  ,   I  tils  is  often 
true  even  when  the  number  of  ships  available  is  limited.  At  first a  it 
would  appear  that  this  is  a  problem  wherein  the  capacity  restrictions  are 
on  the  nodes  of  a  graph  rather  than  its  arcs'.   We  suggests  however s  that 
this  can  be  reduced  to  the  previous  problem  by  employing  a  device  similar 
to  that  used  by  Orden.   See  Chapter  I.  1  sde  representing  a 

port  of  limited  capacity  as  both  a  receiver  and  a  shipper  and  replace  the 
single  node  by  two  nodes  joined  by  an  arc  showing  the  capacity.  We  have 
merely  introduced  a  slight  modification  to  the  model  to  show  more  detail. 

The  applicable  units  of  capacity  might  be  troublesome  in  a  ship- 
ping network  d»e  to  the  wide  variety  of  ships  with  which  one  might  deal. 
There  is  a  great  temptation  to  dismiss  this  difficulty  by  introducing 
some  sort  of  a  "standard"  or  "notional"  ship,  and  for  a  first  approxi- 
;ion  this  may  suffice.  However,  some  j     may  well  be  at      handle. 
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for  example,  3  notional  ships  per  day  provided  this  capacity  actually 
came  in  three  ships  near  our  "s tandari"  s£zes  and  be  incapable  of  even 
docking  one  very  large  ship  of  equivalent  capacity „ 

Thus  our  model  begins  to  break  down  and  we  need  to  Introduce  a  more 
realistic  one  capable  of  handling  the  real  situation,,  We  would  emphasize 
that  the  models  with  which  we  have  dealt  to  this  point  are  mostly  systems 
models  dealing  with  large-scale  inputs  and  gross  effects,,   They  would  re- 
quire inputs  representing  aggregation  of  many  lesser  inputs,  and  these 
must  often  be  combined  by  skillful  application  of  lessons  provided  only 
by  experience o   Simple  arithmetic  summation®  rarely  suffice „   The  inter- 
pretation and  application  of  the  results  provided  by  the  model  is  equal- 
ly difficult,,  and  requires  the  same  careful  Judgment „ 
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MINIMUM  TRANSPORT  UN  •■  A  PD 

lo   General  Discuss ion „ 

We  shall  consider  in  this  chapter  a  problem  wbi'.         a  sense, 
complimentary  to  the  problem  of  maximal  flow  previously  discussed,,  We 
shall  be  concerned  with  the  determinati '     the  minimum  n        trans- 
port units  necessary  to  maintain  a  fixed  schedule o 

In  order  to  make  precise  the  notion  of  a  schedule;,  we  shall  again 
use  the  language  of  graph  theory,,  We  envision  a  connected  graph  in  which 
a  subset  of  nodes  is  distinguished  as  terminals;  others  represent  junc- 
tion pointSo   Sy  a  run  we  mean  the  traversing  of  a  certain  path  from  one 
terminal  to  another  by  a  single  transport  unit.   Runs  originate  i depart) 
and  terminate  (arrive)  at  certain  specified  times  within  an  interval 5, 
such  as  a  day  or  a  week,,   These  departures  and   arrivals  are  repeated  in- 
definitely in  past  and  succeeding  intervals  of  the  same  duration,,   We 
choose  a  typical  time  interval  of  length  77*  9   commencing      7*e  arbi- 
trary time  t   and  ending  at  time  t,  -  t  +  7T  „  A  schedule  is  the 
o  1    o         ■ — ' — ' 

specification  of  the  time  of  each  departure  and  arrival  during  the  time 
interval  Jf   9  for  each  terminal  in  the  system,,   It  will  be  noted  that  the 
duration  of  a  run  may  be  less  than8  equal  : ■•■  or  greats.        a    typical 
interval „ 

Perhaps  the  simplest  feasible  solution  ©f  a  schedule,  problem  is  to 
assign  each  unit  to  repeated  round-trips  betw«     k    terminals,,   This  may 
involve  the  movement  of  a  great  many  empty  or  "dead-head-"  unitss  and  is 
rarely  optimums  if  other  utilization  of  units  is  feasible „   As  can  be 
readily  seen8  the  problem  of  minimizi'  nifl 

scribed  schedule  reduces  to  t  npties  fij  9  C2J  „ 
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We  have  found  several  procedures  for  de  terming  minimum  units  for 
fixed- schedule  operations,,  [,2"]  ,  [.3  J  .   L4J    (jjj         I    *pter  we 
shall  describe,  as  a  representative  example,  an  algorithm  by  To  E.  Bartlett 

£3  J  built  around  a  railway  schedule. 
2,   Algorithm  for  Minimum  Transport  Units. 
Assumptions 

(1)  Schedules  under  consideration  refer  to  runs  of  transport  units 
between  terminals^  The  runs  originate  or  terminate,,  at  specified  times,, 
within  an  interval ,  of  specified  duration,  and  are  repeated  indefinitely 
in  past  and  succeeding  intervals  of  the  same  length, 

(2)  Each  terminal  has  at  least  one  route  connecting  it  with  one  or 
more  other  terminals.  Any  terminal  may  be  reached  from  any  other  in  the 
system.   The  system  may  be  represented  by  a  co     *.d  graph. 

(3)  Each  terminal  has  as  many  "departing1*  runs  as  "incoming"  runs. 
This  may  require  transfer  of  empty  or  "dead-head"  units. 

(4)  For  each  run  there  is  a  standard  or  minimal  type  of  transport 
unit  which  is  employed.   Units  may  be  employed  on  any  run  for  which  speci- 
fications are  met,, 

The  above  assumptions  are  quite  reasonable  and  are  usually  satis- 
fied by  transportation  systems  operating  :  n  .schedules. 
Additional  Assumptions 

For  the  particular  algorithm  presented,,  certain  additional.,  restric- 
tive assumptions  were  made.   These  were. 

(1)  Requirements  at  any  terminal  for  transport  units  to  make  up  out- 
going runs  must  be  obtained  from  regularly  scheduled  incoming  runs  at  that 
terminal. 
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(2)  All  transport  units  are  of  a  single  type,, 

(3)  All  schedules  are  met  precisely,, 

(4)  No  delay  beyond  the  time  differential  of  scheduled  arrivals 
and  departures  is  required  for  interchange  of  transport  units  from  one 
run  to  another o 

Discussion 

By  the  hypothesis  concerning  the  schedule  a  unit  is  either  on  a  run 
or  it  is  "idle"  at  some  terminal  awaiting  transfer  to  a.  run,   A  schedule 
which  is  maintained  generates  a  total  running  time  plus  idle  time  which  is 
equal  to  the  number  of  units  employed  multiplied  by  the  length  of  a  typical 
period.   Since  the  total  running  time  is  fixed  b^  the  schedule^,  the  number 
of  units  employed  will  be  a  minimum  if,  and  only  if,  the  total  idle  time  is 
at  a  minimum.   Therefore,  what  is  required  is  a  pairing  of  arrivals  and 
departures  at  each  terminal  which  will  minimize  total  idle  time, 

Three  separate  cases  are  considered., 

Case  I.   All  departures  in  a  typical  period  at  a  terminal  are 
later  than  all  arrivals  during  this  period. 

Case  II.  Some  departure  is  earlier  than  some  arrivals,  but  all 

arrivals  and  departures  within  the  period  can  be  paired 
in  at  least,  one  way  so  that  each  departure  is  later  than 
its  paired  arrival „ 

Case  III.  Up  to  some  time  during  the  periods  more  departures  than 
arrivals  have  been  scheduled. 

In  the  first  case  any  pairing  of  departures  and  arrivals  generates 
the  same  total,  idle  time.   In  the  second  case,,  pairing  of  the  "first  in" 
with  the  "first  out"  is  both  feasible  and  yields  minimum  idle  time. 

Case  III  presents  a  special  problem.   Since  there  are,  up  to  a  point, 
more  departures  than  arrivals,  there  is  no  feasible  pairingof  all  departures 
with  arrivals  within  one  period.   It  is  necessary  to  pair  one  or  more 
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departures  with  arrivals  I    a  prece      .      and  seme  arrival  (s) 
within  the  period  are  lei 'f  jvet  for  pairi  -      depart     lin  a  sub- 
sequent period.   By  extending  the  sequence  I  ith        rward  and  backward 
sufficiently  to  pair  all!  departures  wi  ig  problem 

reduces  to  that  of  Case  II,  had    ■  Lrst-in-f      at  proc      gives 

minimum  idle  time.   However,  it  is  r- :  -         idcpt  a  consistent  assoc- 
iation of  idle  time  to  a  particular  period  in     i  not  to  miscount  the 
total  time.   The  procedure  adopted  was  to  associate  with  the  current  per- 
iod that  idle  time  generated  by  pairing  nans  of  the  previous  period  to 
runs  in  the  current  period  as  well  as  the  idle  time  generated  by  pair- 
ing within  the  period,, 

Idle  time  developed  in  the  pairings  from  the  current  period  to  the 
next  period  is  not  charged  to  the  current  period^  but  to  the  next  one„ 
Results 

From  the  above  considerations  a  tot  1   i  lining  plus  idle,  time  is 
calculated,  and  is  shown  to  be. 

T  -  It  (  0(  +  %  fi  . ) 

»re 

if   s  length        -d 

oC     =  total  runs  I  lining  (or 

end)  of  a  pea     !  i  all  tei  ninals 

tD'    -   maximum  value  of  cumu]    we,   departures  less 
cumulative  arrivals  at  terminal  j  during  a 
period,, 

The  minimum  number  of  transport  units  required  is  then  obtained  by 

dividing  T  by  the  length  of  the  period  *Jf  ,  and  if  given  as 

ID  =  minimum  number  of  transport  units  ~  c(   +  j>  6 
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In  concluding  his  paper.,  Bartlett  pointed  out  areas  of  schedul- 
ing problems  under  study.  I  can 
be  seen,  are  aimed  at  removing  the  restrict  s  made  earlier. 

Items  listed  were; 

a.   Determination  of  routings  to  permit  economic  maintenance 
schedules. 

b„   Determination  of  routings  to  optimize  the  number  of  units 

used  for  the  situation  in  which  several  classes  of  equipment 
may  be  employed  on  various  routes. 

c.  Determination  of  routings  to  minimize  the  number  of  utti 
used  and  maintain  schedules  having  probabilistic  arrival 
t  ime  s . 

d.  Determination  of  optimum  schedules  based  on  cost  of  trans- 
port units,,  operating  costss  and  demand  distributions. 

e.  Analysis  of  the  possible  application  of  the  concepts  : 
the  field  of  machine  scheduling. 

The  actual  problem  of  allocating  units  to  a  schedule  in  a  manner  to 
preserve  minimality  ands  at  the  same  time?  meet  certain  maintenance  require- 
ments, was  discussed  in  a  su  bv  T.  E.  Bartlett  and  a.  Charnes 

p>J  ,   Their  method  of  solution  involves  translating  the  schedule  into  an 
oriented  or  directed  graph8  and  working  wiJ  Lated  incidence 

matrix.  A  node  in  the  graph  considered  represents  a  specific  time  it  a 
specific  terminal.   The  time  chosen  is  the  time  of  a  local  maximum  in  the 
"idle  equipment  inventory"  at  the  terminal  under  consideration.   There  may 
be  more  than  one  "node"  for  each  terminal;  there  will  be  at  least  one.  An 
arriving  run  prior  to  a  node  time9  which  may  be  paired  with  a  departing  run 
after  the  node  time-.,  without  violating  minimality  of  equipments,  is  regarded 
as  a  link  (arc')  positively  incident  on  that  node.  A  departing  run  after 
a  node  time  which  can  be  paired  with  an  arriving  run  prior  to  the  node 
time  is  viewed  as  a  link  negatively  i  jn  the  i  ■ presenta- 

tion gives  sn  .   Idea  cf  the  versatility  of  gj     :heory  in  the  study  of 
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transportation  problems 

With  the  incidence  matrix  constri  ".  ribed  above,  j,  each 

column  corresponds  to  one  of  the  runs  of  the  system9  and  each  row  corres- 
ponds to  one  of  the  nodes „   One  can  then  trace  out  a  cycle,  or  a  "pool", 
in  the  following  manner,,   Connect  a  minus  one  with  a  plus  one  in  the 
same  rows  and  repeat;,  connecting  plus  ones  and  minus  ones  alternately 
within  the  same  column  and  then  the  same  row0   When  the  first  or  start- 
ing entry  is  connected  the  second  time,  a  cycle  is  formed,   The  pairing 
of  runs  in  this  manner  is  both  feasible  and  maintains  minimality  of 
equipment „ 
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CHAPTER  V 
MISCELLANEOUS  TECHNIQUES 

Up  to  this  point  we  have  dealt  with  a  rather  restricted  group  of 
transportation  problems  and  models.   For  the  most  part  the  models  have 
been  of  the  network  type  and  descriptive  of  an  entire  transportation 
system   In  this  chapter  we  shall  give  examples  of  the  application  of 
other  techniques  to  problems  in  transportation.   Specifically  we  sh 
cite  examples  using  queuing  theory,  the  Monte-Carlo  method,  and  engineer- 
ing analysis. 

Among  the  earliest  models  employed  in  the  study  of  transportation 
systems  was  the  queue,  or  waiting-line,  model.  Examples  of  queues  in 
transportation  systems  might  include  automobiles  at  a  stop  sign  or  a  toll 
station,  freight  cars  in  a  freight  yard,  airliners  circling  an  airport 
awaiting  clearance  to  land,  or  ships  in  harbor  awaiting  unloading  In 
a  different  but  related  field  this  model  has  been  applied  to  telephone 
calls  at  switchboards  and  messages  in  communications  relay  stations. 

Waiting  line  theory  is  strongly  related  to  probability  theory, 
since  both  the  arrival  times  of  "customers"  and  the  service  time  for  an 
individual  customer  are  stochastic  variables   This  is  true  to  a  large 
extent  even  though  arrivals  might  be  "scheduled".   The  usual  practice  is 
to  study  the  so-called  steady  state;  that  is,  one  assumes  the  queue  has 
been  in  operation  long  enough  for  transient  effects  to  have  diminished  to 
a  negligible  point ;   One  usually  knows,  or  can  soon  estimate  with  accep- 
table accuracy,  the  mean  arrival  rate  of  customers  and  the  mean  service 
time   The  specification  of  the  distribution  of  these  random  variables 
is  more  difficult   However,  the  most  widely  used  assumption  is  that  the 
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number  of  arrivals  within  a  specified  interval  has  a  Pol     distributions 
and  that  service  times  have  a  negative  exponential  distribution,  Wi 
these  assumptions  and  inputs,,  the  model  can  be  analyzed  to  predict  mean 
holding  time  -  waiting  time  plus  service  time  -,  mean  quev  !      t9  and 
probability  a  customer  will  have  to  wait, 

A  good  example  of  the  application  of  queuing  theory  to  a  transport 
tion  problem  is  contained  in  an  article  by  Leslie  C,  Edie  Qj  of  The  Port 
of  New  York  Authority „  in  which  he  discussed  traffic  delays  at  roll  booths 
of  several  tunnels  serving  the  city  of  New  York,, 

In  ■     problems  involving  queues,,  solutions  in  closed  form  are 
difficulty  if  not  impossible,  to  obtain,,   Such  problems  are  sometimes  amen- 
able     certain  amount  of  analysis  by  the  Monte-Carlo  method.   In  using 
this  method^  one  develops  a  model  of  the  system  he  wishes  to  study  and 
then  simulates  arrivals  from  time  to  time.   The  arrival  times  are  assign- 
random  in  accordance  with  previously  observed  arrival  statistics, 

model  the  investigator  can  determine  the  effect  of  varying  certain 
parameters  in  the  system.   He  cans  for  instance ^  introduce  an  additional 
service  channel  or  allow  priorities.   With  the  simulation  technique,  the 
researcher  can   get  data  corresponding  to  months  of  actual  operation  quick- 
ly and  cheaply, 

Roger  R0  Crane?  Frank  B.  Browns  and  R„  0,  Blanchard  C2J  have  reported 
on  the  application  of  this  method  to  a  railroad  classification  yard.   These 
authors  developed  a  model  by  the  method  described  above  ands  using  simpli- 
fied statistical  assumptions,  isolated  areas  for  more  careful  study.   In 
the  two  areas  reported  they  determined  the  effect  of  substituting  parallel 
processing  for  series  processing  and  of  changing  the  "rules"  regarding  the 
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nes,   Considerable  saving*  were  realized  in  both  in- 

in  this  context  c      ed  of  insp     i  and  cla 
fication. 

For  our  last  exam;,      nave  chosen  the  Cargo -Hand  ling 
Froject  in  the.  Department  of  Engineerings  University  of  California  at 
Los  Angeles,  This  research  is  sponsored  by  the  Office  of  Naval  Research 
and  the  Maritime  Administration,   Though  by  no  means  limited  there;: 

:  project  has  made  use  of  engineering  analysis,  statistical  methods,, 
and      onte-Carlo  technique.   The  project  has  been  reported  in  detail  in 
a  series  of  numbered  reports^  (references  14  J  ,  L5 J  9    \_6j    ,    [J J    i,   rhese 
arc;  not  generally  available,   R,  R,  0JNeill  reported  on  the  project  a' 
discussed  the  analysis  and  Monte-Carlo  simulation  of  cargo  handling  be- 
fore the  Purdue  University  First  Transportation  Research  Symposia 
February  1957,  This  presentation,,  along  with  others  given  at  the  sympo- 
..  was  published  in  Naval  Research  Logistics  Quarterly  £"3j     -pt ember 


s  transfer  of  cargo  between  ship  and  shore  has  long  been  recognized 
as  the  major  bottleneck  in  the  operations  of  the  shipping  industry.   The 
basic  objec.-i.-ct"  of  the  research  at  the  University  of  California  at  Los 
Angeles  has  been  the  effecting  of  improvements  to  the  cargo-handling 
system. 

For  the  analysis  of  cargo  handling,,  the  project  faced  the  problem 
of  choosing  a  level  of  description  which  would  yield  to  analysis  and  yet 
lead  to  useful  results.   The  level  chosen  considered  the  cargo-handling 
system  as  a  series  of  transporting  links  in  which  each  link  was  the  sum 
of  all  the  movements  required  of  a  carrier  in  transferring  a  load  of 
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cargo  from  one  point  to  another,,   In  a  land-to-water  sequence  the  links 

studied  were:  {a)      the  shed  link,   (b)  the  wharf  link 

link»   (d)   the  hold  intermediate  link,,  and   (e)   the  hold  link„ 

With  the  cargo-handling  system  idealized  in  the  foregoing  manner 

certain  hypotheses  were  made  concerning  the  interrelations  existing 

among  the  various  links.   Two  examples  of  such  hypotheses  are: 

Hypothesis  I  If  there  is  a  group  of  transporting  links  connected 
in  series,  then  there  is  one  link  which  has  the  least  mean  induced 
delay  and  is  the  slowest  or  controlling  linko 

Hypothesis  II:  The  relationship  between  the  relevant  factors  which 
affect  the  mean  activity  time  of  a  link  can  be  formulated  so  as  to 
make  prediction  possible,, 

Both  of  these  hypotheses  were  used  as  the  basis  for  a  field  study  con- 
ducted in  the  summer  of  1953s  primarily  in  the  ports  of  Los  Angeles  and  San 
Francisco „   Even  though  a  randoms  work-sampling  technique  was  used  to  obtain 
data,  there  were  41,218  individual  observations  recorded*   This  gives  some 
idea  of  the  magnitude  of  data-gathering  problems „ 

It  is  noted  that  at  the  gross  level  a  port  terminal  can  be  thought 
of  as  a  node  with  a  land  carrier  as  the  link  on  one  side  and  a  ship  as 
the  link  on  the  other  side0   The  approach  used  in  cargo-handling  study 
Illustrates  the  meaning  >f  a  remark  made  earlier  concerning  the  necessity 
of  buildinj.  models  in  a  degree  of  detail  consistent  with  the  problem  to  be 
studied,, 
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APPENDIX  A 
ALGORITHMS  FOR  THE  HITCHCOCK  DISTRIBUTION  PROBLEM 
1,   The  Dantzig  Algorithm, 

In  this  appandlxj,  we  shall  describe  the  essential  elements  of  Dant- 
zig" s  algorithm  for  solving  an  m  x  n  transportation  problem  based  on  his 
simplex  method,,   See  Chapter  I,   We  shall  then  work  out  a  simple  problem 
using  this  method.   We  shall  make  similar  presentations  of  the  algorithms 
of  Ford  and  Fulkersons  and  Munkres,  working  out  the  same  problem  in  each 
instance,,   The  problem  is  4  x  5  (four  sources,  five  destinations)  and  the 
numbers  are  quite  simple.   We  hoped  to  choose  one  which  could  be  followed 
easily  and  yet  was  sufficiently  difficult  to  illustrate  the  methods,, 
problem  is  stated  below° 

There  are  four  sources  with  amounts  to  be  shipped  as  follows" 

•x-w 

a2.60 

a  =  40 
a4  =  90 

There  are  five  destinations  with  amounts  required  as  follows; 


bl 

- 

20 

b2 

= 

40 

b3 

= 

70 

\ 

= 

30 

b5 

- 

40 

We  note  5   a.  =  5  b.  =  200. 
i        J 
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The  cost  of  transporting  a  unit  amount  of  commodity  from  source  i  to 

destination  j  is  given  as  the  ij  entry  in  the  following  matrix: 

_1    Z  3  4- 

/ 

Z 

3 


4 

l 

6 

5 

/ 

Z 

3 

jT 

3 

Z 

5 

z 

3 

4 

* 

4 

1  \  * 

1 *— ■ — — 

1 

2 

Our  problem  is  to  findc 

X  =   Tx„  ,~l  such  that 


x  .  >  0 


x  ,  -   a , 


5"  x. .  =  b 

2fc 


<•  3 


.  x,  . 


i  =  1,2,3,4  j  =  1,2,3,4,5      1.1 

i  -  1,2,3,4  1,2 

j  -  1,2,3,4,5  1,3 

=  Minimum 

Before  commencing  the  algorithm,  we  shall  state,,  without  proofs  some 

useful  theorems  relating  to  the  probiem„   For  proofs  the  reader  is  referred 

to  Dantzig  CO  and  Gass  C^2   . 

Theorem  1„   The  transportation  problem  has  a  feasible  solution. 

Theorem  2.   A  solution  of  at  most  m  +  n  -  1-positive  x,  9s 
exists.  J 

Theorem  3.   If  the  a,   and  b   are  all  non-negative  integers^ 

then  every  basic  feasible  solution  has  integral  values, 

Theorem  4.   A  finite  feasible  solution  always  exists. 
For  a  more  general  case  with  m  sources  and  n  destinations,  if  we  write 
out  the  equations  1.2  and  1.3  we  have  m  +  n  equations  in  mn  unknowns.   But, 
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.-.ecified   Ta  =  S   h.  rndant  equati 

and     system  may  be  reduced  to  m  +  n  -      learly  indpenden: 

A  basic  feasible  solution  then  has  m  +  n  -  1  positive 

feasible  solution  is  found  with  less  than  m  3- n  -  1  positive   x^  °  i 

s  indicated.   The  difficulty  can  be  averted  by  perturbing  the 

following  manner^   add  a  small  positive  amount  £"  to 

leave  all  b.'s  unaltered  except  for  b  9  and  add  n£  to  b  , 
j  n  n 

s  proce     need  not  be  employed  until  degeneracy  is  indicate 
lem  for  a  minimum,  basic  feasible  solution,  using  th«  & 

r   the  solution  is  obtained  allow  each  £  to  go  to 

r.lido   In  our  problem9  since  each  a^  and  b.  is  a  pc 

i      j 

tive  integer,,  we  are  assured  of  a  solution  involving  in- 

tegers. 

To  explain  the  computational  procedure  we  shall  work  out  our  p  ; 

the  initial  step9  we  write  down  the  matrix  of  direct 
cost  coefficients  and  a  blank  assignment  matrix.   We  shall  habit uai 

Lumn  of  a.'s  to  the  right  of  the  assignment  matrix,,  and  a  row 

of  b  *s  below  the  matrix,.  This  procedure  will  serve  to  keep  the  pi 
j 

m  before  us» 
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5 

2 


1 


L      S-     I 


- 


5 


4 


I 


dc 


10 


' 


b:  i  Zo   fe    7o   zo  4-o  ' 

J       ! -4 J I  • 


th  the  as-      . 
the  minimum     f_  '  .J  ]  0   ui  - 

jrther  cell    (231)   and 


b,    -    10)    =   10  We    have  now 

1  ' 


complete      ignments  in  column  onet   We  next  move  to  cell       and 
nl  minimum  L40j,60  -  10  J   =  40„   Continuing  in  t 
In  the  following  basic5  feasible  solution,, 


/ 


J 


<\ 


/ 

L. 

•> 

4 

5 

lo 

lo 

4-0 

*ol 
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In  this  example  we  have     m  +  n  -    1   =  8   p      Ltive     x 


;  f  aft 


having  made  one  of  our  assignments  x       amount  reri 

umn  j  and  that  remaining  in  the  row  i  had  /^rushed  simulta--  ^e 

Id  have  had  a  degenerate  solution,,  except  for  the  case  i  -  d  j  = 

Bo   In  case  of  degeneracy 8  inclusion  of  the  £*  perturbation,  as  described 

above,  will  remove  it;  for  then  no  partial  sum  of  the  a  u s  and  b  „ " s  can 

equal  any  a.  or  b,. 
'     J 

For  comparison  purposes  we  note  here  the  cost  of  our  first  basic 

feasible  solution,, 

C,  a  10(4)  +  10(2)  +  40(3)  +  10(5)  +  40(3)  +  20(4)  +  30(1)  +  40(2)  ^  540 
1 

-xt  we  computfij  for  the  basic  feasible  solution  just  obtained,  its 

indirect  cost  matrix,,   Each  element   c,  .  of  the  indirect  cost  matrix  will 

ij 

be  the  sum  of  an  indirect  cost  u,  associated  with  the  source  (row)  and 

an  indirect  cost  v.  associated  with  the  destination  (column).   Thus 

J 

for  all  i  and  j„ 


ij 


u.  +  v . 


For  those  pairs  i .,  j  employed  in  tht  trial  solution 


c„  .     c.. 
ij      ij 

From  thus  latter  equation  we  can  get  some  of  the  entries  for  the  indirect 
-  cost  matrix.   In  our  example  we  get  the  incomplete  matrix  below. 

rect  -  Cost  Matrix  (incomplete  i 


1       —  ■ 

4  ! 

2 

3 

5* 

3 

4 

/ 

Z 
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We  ••.. 


nd  v.  =  Os  art 


J 


for  a] 
Lndirect  cost  mat 

:e  by  the  eq        uu  4  v.  -  c_  -  2^   u  •  2„   Then  v^  =  1 

and        until  all  uu   and  v0  at'  are 

1        j  ij 

relations   c. „  =  u.  + 


Indirect  Cost  1 


vi 


U; 


4 

1 

7 

S 

2 

5" 

3 

4: 
/2 

z 

3 

0 

/ 

3 

0 

/       0 

0 

2 

- 

t 

z 
1 

1 

f   > 

3 

o 

f 

i 

order  to  fat     ce  comparison  :nding  entries  in  the 

ditect  and  indirect  cost  matrices*,  we  shall  write  them  together  as  one 

called  the  cost  matrix,,   In  each  cell  of  the  cost  matrix  we 
sha  in  t  and  c,  ,  in  the       ialf„   In  ■ 

sample  problem  this  res  ii owing  cost  matrix  for  seep  i 
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vi 


Cost  Matrix 
(Step 


Ut 


4  ' 

9[ 

5   \ 

n .  ■ 

'■  (o 

/  2- 

4 
!  2 

2 

2. 

5 

3 

2 
,'3 

0/ 

1 

Z 

5' 

0-' 
/4 

/4 

0 

/ 
4 

z  ■ 

/4 

r       / 

0 

I 

3 

r 

o 

/ 

Entri*:        the  diagonal  are  c.  .    L  =  1 

1J   j  -  1 


aal  are   c 

it 


■  .  \ 


*■  J 


j  =  1 

]    i     •  \J   - 

The  additional  row  I        ^.he  row  of 

included  for  conver. I 


I  '' 


We  now  look  at  those  cells  in  which  c   >  c   „   If  there  are  nones 

the  basic    ^slble  solution  we  have  chosen  is  a  mim~  and 

our  problem  is  solved,,   If  c   >  c  ..   for  some  i  and  ja  w< 

determine 

M  =  max  c ,  .  -  c 

ii  lj      1J 

>ur  case  M  =  49  for  cells  (ls2)  and  1*4)     ire  ost 

step  I,   We  wish  now  to  assign  to  one  cell  corresponding  to  Ms 

1 
--•  -.-■.  0        as  large  as  possible  consistent  with  the  requirement 


.1  x.  ,  be  ncn- negative.-,  and  that  S    x. .      -     a.       2E     -  b,. 
If  we  increase  the  assignment  in  cell   (i  ,  j  )  by  v.9   we  must  decrease- 


by  a  corresponding  amount  the  assignment  in  some  other  cell  of  row  £ 
order  to  meet  the  Hatter  constraints  above.   This  adjustment  will  cause  an 
unbalance  in  a  different  column;  so  we  must  continue  this  compensating 
procedu     .til  we  complete  a  cycle  by  decreasing  by   ©  the  assigi 
in  some  cell   (i,  j  ). 

e  two  candidates  for  the  assignment  ©  „  cells 
We  arbitrarily  choose  eel]     I)   and  formally  enter  0 
■i  -,  vi  lu    Making  the  required  ad- 

jus*-  nt  matrix  in  the  following  interim  form. 

cycle  we  have  foil      alternately       and  subtracting  © 

J. 

I 


1   The  subscript  I  refers  to  the  first  iteration.   In  subsequent 
iterations  we  use  O^     O    „  . .  ,  until  the  problem  i.; 
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1 0-6, 

-, 

/ 

40-0, 

10 

<?o 

2o 

3d 

*\Q       lo 

b-  i 

2o 

- 

10 

3o 

In  each  iter.'      at  this  poit        will  alway 
entri<     f  x.   -  0  J    We  make  ©  as  large  as  pos 


1J 


x        >  he  present   example  0 ,     ^      10  ii  x, , 

n-ne      tive,         i .     ■  •:    set       @  .   =   10  arid        .  net  rix 


Assignment  Matrix 


[step   2: 


U.L 


k 


1 


l 


/o 


30      /<? 


r 


/6 


ZO  |40, 


» 
7>)     3^!4^ 


Lted   In  a  new,    basic   feasible   solution 


p     iiti  id     x  -  Oo      If 9    at   any   stage s    tl 

to  one   route   causes   two  or  more  positive     x,         to  varus... 

ij 

'■ndicated.      The   inclusion  of   the     €  perturbation,,    pre 

ds    will   eliminate  the  degeneracy, 

47 


Low  we  rewrite  our  assignment  matrix  (step  2)  and  i       elated 

Step  2 

Assignment  ; 


Cost  Matrix 


U, 


0  -' 
'4 


/ 


4 


M    * 


/ 

,  2 


j 

5 

i 


si/ 


4,' 


C       Z 


'4 


I 


/ 


/ 


■ 


i    - 

4 


/-? 


1 


-  /  f  C 


a. 

^ 

/0 

4o 

*> 

0, 

.      , 

5 

4o 

hiave  again  encircled  two  possible,  routes  for  the  introduce.      0- 

Again  car  choice  is  somewhat  arbitrary,  but  route  (4f2)  has  t1 

direct:  cost;  hence9  we  choose  it„   We  have  indicated  the  plus  and  minus 

(3)o's  above.   From  the  assignment  matrix  we  see  0~  ^  20^  in  order 

that  x,   not  be  non-negative.   We  let  @  „  «  20,   Continuing  this 

process,  we  go  through  the  following  iterations0 

Step   3 

0.1 


"-"•l 

0 

1 

/ 
1 

3 

o 

/ 

>2 

/ 

/ 

Z 

J 
3 

O 

3 

■ 

0 

b 

O 

4 

/ 
/  Z 

I  ■ 

3,    ' 

/4 

'  4 

i 

1 

rrA 

J  I 

-1 

t 

0 

1 1 

2 

4 



1 

Zo 


10  \ 

I 

j  4o 


30 


/0 

03     bo 

9o 


(%,Jo 
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Step     4 


tL\ 


U^i 


V\ 


2 

4 

• 

76 

■  —  ■-> 

/ 

: — i 

z 

\ 

y  . 
/2 

I  ■ 
.'3 

S~ 

/3 

z/ 
/'z 

/ 

0 

-I,' 

-1  / 

0/ 

-/ 

2 

■ 
l   ■ 

s 

■  4' 

1  , 
i 

/z 

/ 

1 

0 

4 
1 

1 

1 

10-ty 

Oo 

! 
/o' 

Zo 

So 

10 

t',0 

4o 

! 

3o*% 

3o 

36-3 

9o 

Zo 

40 

70 

3o 

- 

Step      3 


UL, 


V_ 


/  1 

■    - 

6 

-, 

r    -        —, 
1      ' 

/      0 

i 

z  -' 

'Z 

/  - 

,  '5 

1/ 

/  3 

2/' 

■  i 

z 

0 

-t  7 

/  z 

3  7 

/3 

- 1  ' 

o/ 
<  4 

Z/ 

1  / 

'  1 

••  V  ■ 

1 

o 

-/ 

3 

1 

- 1 

0 

1 

0<t*io 


Ui 

10 

10 

2o 

3o-e>i 

/o>0s 

\(o0 

4o 

lo 

0s- 

3o 

10- Qs 

9o 

4 

Zo 

40 

70 

i , 

10 

0<r  ^Zo 


Step     6 


u: 


vl 


/ 

4/ 

-6 

■  i 

i  / 

1 

/Z 

/3 

> 

Z'' 

,'z 

2 

o/ 
-'  s 

• 

t>  ■ 

• 

■S  , 

o/ 

0 

1 

A 

1  / 
/ 1 

/4 

/    ■ 

t 

1  / 
/z 

1 

_/ 

0 

3 

0 

0 

LX-t 


' 

10 

to 

\Zo 

10 

30 

10 

io 

Zo 

3o  \ 

L 

9o 

— __< 

Zo 

4o 

70 

30  ■  ,o 

, 
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-    step  6,-,   M  s  max.  (c,  .  -  c.  .)  =0,  and     issignment  ma" 

i  J 

represents  as         I he  cost  : 


i  ■. 


[  3 )  +  4 


+  30 1 

»  10  *  40  +  50  +  60  +  120  +  40  +  80  4 

2„   Tl         I  tkei    Algorithm, 

Two  matrices  are  carried  along- 

(1)  The  cost-weight  matrix.   This  is  the 

with  "\  sched  to  the  rows  and 

An  additional  matrix  of  the  same  size  as  C9  call 

matrix,   rhis  matr      «sed  in  making  the  assignment ss  i0e0 

the  positive  x  's.   With  each  row  we  associate  a  "surplus"  and  wi 
■*■  J 

each  column  a  "shorta--     --ttially  these  are  the  a  ani  b  of  the  original 
problem,      Ldez  the  problem  previously  given  whi<  i   in  the 

following  tableau. 


4 

/ 

b 

5 

/ 

2 

3 

S~ 

3 

2 

S 

Z 

3 

4 

4 
2 

4  !    / 

4 

/ 

Zo 

4o 

70 

SO  I  ^o 

10 

6o 

40 

9o 



>.   obta'  il  cost-weight  matrix  and  :    i  litial  zero 

i  c   Prepare  two  matr 

matrix  C  =   [  c ,  ].   ; 


ij 
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?pt  tc  -.rial  surpluses  a. 


u                        i 

1 

5' 

/ 

\   1 

'■ . 

70 

■ 

, - 

_) 

-1 

4 

/  6 

3 

O 

\-z 

- 

3 

. 

1 

/>      5 

4 

1     4 

6 -  -M 

1 — U- 

A 

! 

9o\ 

o 

Q 

of  the  cost  matrix,  pi  ; 
its  I  >.e  weight  for  that  row,   (ThJ 

in  the  zero  matrix 
the-     st  values  in  the  cose  matrix,, 
abo-:  - 

r'  column  in  the  z 

ass  E  zero  and  enter  these  weights   in  ippi 

a  ro*/  above  the  cost  matrix.   For  each  -  a 

ee  in  tl     imple,  sea  tnatris  b  .  . 
;!  it   .Take  the 
ti  column  weight  for  t  I  to 

.    ponding 
1 1  h  row  and  each  column 
■  m  I    ''tt   matrix*   It  wil] 

-■.own  be  "  ■  I  e  „ 
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-  / 

0 

o 

-/ 

i 

4 

( 

6 

S     1 

-I 

z 

3 

3 

1 

3    Z 

-I 

s 

z 

4     4 

-/ 

1 

4 

3 

\ 

2.0 

'-{0 

70 

3c 

10 

o 

o 

\6o 

o 

* 

10 

OQ^ 

'■ 

do 

D 

o 

nay  be  as  partially  satisfyinj 

in      i'lumms  from  surpluses  in  the  rows.,  using  only  circled  posit  i 
as  ■;        shipping  routes.   This  step  is  accomplished  as  follows: 

first  row  for  the  first  circle  appearing  in  it.   Determine  the 
he  surplus  for  that  row  and  the  shore      n         imn. 
Enter  this  number  in  the  circle  and  simul1    »usly  reduce  the  shortage 
surplus  by  thai    i  it0  Then  i  cond  circle  and 

'/hen  the  first  row  is  completed       the  second.   Repeat  until 

been  scanned.   Applying  this  procedure  to  the  zero  matrix  of  our 
problem  yields i 


\ 

o 

o 

' 1 

6o 

o 

o 

■ 

\o\ 

® 

o 

© 

o® 

I  ^^ 

.- 

Go 


:     stive    ; 

lit 
or  non-br 
discussed    >e] 
label: 

by   iabeling   the  mple 

tl  e    label      i  JL  £®    JL 

rplus   in  that    row,,  ir- 

■ 

sxample 


\l   o 

'. 

© 

® 

Of) 

\£) 

(Jo) 

O 

. 

© 

j 
J 

1    /* 


ereafter  t  i  C  will    -  sntial    . 

ass< 

.    ; 

point,   We  take  each  ...  circle 

in  currently  unlabeled  columns „   If  .-  i   label  the 

coluira      Lning  it  where  f  is  flow  a     ited 

with  the  rev  being  scanned  and  R.   idem     s  the  row,   Results  of  thi 

x 

:m   be low o 


6-v"      u>Ki> 


C&O.C   ) 


The  number  £  will  hereafter  be  called  the        al  flow'*  assoc- 
iated with  the.  column. 

When,  the  rows  have  all  been  scanned ,  turn  to  the  newly  labeled  columns. 

Scan  each  labeled  column  for  circles  contai  entries  but  not 

appearing  in  a  previously  labeled  row„  *  Any     Ln  which  such  a  circle  lies 

is  given  a  label  {JL     9  C  ,)  where  C,   is  the  column  being  scanned  and 

J         J 

.JL    is  the  smaller  of  the  potential  flow  of  the  column  and  the  entry  in  the 
circle   Column  scanning  differs  from  rou  waysc   The  circ 

x   have  positive  entries,  and  JL    is  computed  different ly„   Applying  this 
step  to  our  example,,  yields  the  following  diagram.! 


(Mu)       ?,**) 
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Having  scanned  all  newly  labeled  columns  we  tarn  again  to  the  rows 
scanning  the  newly  labeled  ones  as  before. 

scanning  newly  labeled  rows  and      as  unci!  either  a  breakthrough 
achieved  or  no  new  labeling  is  |  scant;      u  #3  we  achiev- 

ed a  breakthrough  (i,e,  we  labeled  a  column  that  had  a.   shortage).   At  chat 
point  the  matrix  was  labeled  as  belw. 


.  o    ) 


{^0,00^,^10,%) 


Breakthrough 

We  stop  labeling  when  a  breakthrough  is  achieved  and  make  adjust- 
ments in  the  allocations.   Suppose  we  have  labeled  a  column  which  has  a 
shortage,  s,  with  the  label  (f9  R0)„   In  our  case  s  s  60  and  the  label 
is  (309R  ),   This  means  that  we  can  increase  the  flow  to  the  labeled 
column  by  an  amount  h  =  min  [s»fj  »   The    ':.«:»wing  procedure  is  used, 

(1)   Decrease  the  shortage  by  an  mount  h,   (30  in  our  case) 

Look  at.  the  label  (fsR,  I  a-  i  increase  the  entry  in  the  circle 

X 

where  row  i  and  the  labeled  column  intersect  by  h  units, 

(3)  The  label  on  row  i  designates  a  certain  column,  j„   Decrease 

the  entry  in  the  circle  where  row  i  meets  column  j  by  h  units, 


55 


(4)  Obtain  a  new  row  from  the  label  on  column  j  and  increase;, 
continue  alternately  increasing  and  dec  i  sents  in 

circles  by  h  units  until  a  row  is  i-  ■    label 

■Co>- 

(5)  When  the  indicated  surplus  is  reduced  by  h  units,  this  process 
is  completed c 

Applied  to  our  example  these  steps  produce  the  following  matrix. 
Arrows  indicate  the  "path"  we  followed  in  the  adjustment  process. 


\ 

c 

■ 

r 

--30 

0 

fs 

1© 

i 

& , 

1 

© 

30iio\ 

! 

1 

loO-iO 

\ 

^ 

*50 

L    .    . 

„    IO 

1      1 

If9  at  this  pointf,  all  shortages  and  surpluses  are  zero» 
problem  is  solved.   Otherwise.,  remove  all  previously  applied  labels  and 
commence  the  labeling  process  anew. 

Next  time  around,  we  exhaust  all  possible  labelings  without  achiev- 
ing further  breakthroughs  and  have  the  matrix  shown  below. 
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0 

o 

% 

1 

30      0 

0 

•to) 

0 

0 

20  ) 

0 

s    J 

s 

30 

5o) 

Cz   ) 

(lo ,  d  ) 


'/6,JPt)(5»;jei        ft  ,/'v^ 


Although  we  have  labeled  all  except  one  we 

have  not  achieved  breakthrough,,   At  this  point  we  t  .  .on 

to  the  cost-weight  matrix  which  provide  s  a  means  of  proceedir 

Non-breakthrough 

We  now  change  the  weights  on  the  cost-w«  matrix  by  subtracting 

k  units  from  the  weights  of  the  label)       and  adding  k  units  to  the 

weights  of  the  labeled  columns,,   The  value  of  k  is  determined  by  taking 

it  as  large  as  possible  subject  to  the  constraint  c „ .  +  w,   +  w  >  0 
6     v  J  i     ] 


where  c    ;s  the  ij  entry  in  the.  cost  matrj 


ij 


w„   is  the  row  weic 
i 


w.   is  the  column  weights  The  determination  of  k  can  be  accomplished 
as  follows 

Take  up  the  cost-weight  matrix  but  ignore  for  the  time  being  all 
labeled  columns,.   Shading,  covering,  or  some       list mguishing  mark- 
ing is  recommended  for  the  labeled  columns,,  For  •        y  in  each  label- 
ed row,  compute  the  algebraic  sum  of  it,  its  row  weight,  and  its  column 
weight .   The  minimum  value  obtained  this  way  is  k,,   Enter  circles  in  the 


o  matr  <ositions  corresponding  to  t  apply  I  tep 

s   denote    lab< 


Mat r ix 


c 

A 

/ 

-  / 

-/ 

5" 

l 

-.2 

5 

^2       3 

Z 

... 
-2     f     2 

1 

r3      ^ 

4 

-1 

4 

1 

..S 

1 

/2    ' 

.   Matrix 


o 


o 


o 


iZo 


o 


5G 


o 
fo 


30 


3o 


o 


o 


O 


4  ) 


2 

ase  all  eo    •'  have  '  .We  need  con- 

sider i        tions  .        ,3)  j  and  «,  •      '•■  re  placed  in 

positions     I  and  (4,3)^  indicated  temporarily  by  diamonds  in  the  zei 

matrix.   Subtract  k  from  each  labeled  row's  weight.,  and  add  k  to  each 
labe  ' s  weight.   Now  looking  at  the  .     latrix;  we  ignore  all 

rows  and  remove  all  circles  remaining  in  lal  :  ily 

i  in  pc       (3s, 2)  is  our  cas<  cross- 

:ve,     a  check  on  the  proc-  ,  id  ne'> 

lead  to  the  t  i  circle  with  a   \  '©.  entry., 

above  procedures  provide  us  a  new  zero  matrix  with  which  we  re- 
commence the  labeling  process „   A  new  cost- weight         i   also  shown „ 
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-3 

4 

4< 

i 

■ 

/   U* 

5        1 

-4 

2 

3    5"» 

3     2 

-2:|  5~~ 

2     3 

t  4 

4 

~3  1    4 

i 

/    i^2 

1 

■ — 

With  very  little  labeling  we  •  ough  into  column  3, 


\ 

0 

o 

3a 

0 

o 

0 

© 

o 

o 

£) 

0 

© 

J 

© 

3o 

t— 

© 

o 

® 

(30,  Co) 


M)l 


3o;%;  (54  ^ 


This  breakthrough  provides  us  the  needed  30  units  flow  to  column  3,  and 

the  problem  is  solved*   All  shortages  and   surpluses  are  now  zero* 
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.-  ion  is  given 
■  ■ 

C  =  430 
This  is  the  saint  i ned  by  the 

The  above  description  is  really  more  Involved  than  the  act  ■    se  of 

the  algorithm^   Once  the  ground  rules  are  understood,  the  manipu  I   i  of 

the  matrices  is  not  difficult   The  authors  of  thi 

completed  a  20  x  20  optimal  assignment  problem  in  about  30  minutes 

hand  computation*   The  same  problem  by  the  simplex  method^,  required  well 

over  an  hi. 

3c   The  Munkres  Algorithm, 

Because  of  its  basic  similarity  to  the  For'!  -  Fulkerson  method,, 

s  algorithm  will  be  stated  completely,  and  :     the  example  will  be 

worked  out.   Munkres  statement  of  the  problem       an  to  speficy  hus 

notation. 

Let-. 

D  =  [d„ J  be  an  m  x  n  matrix  of  non-negative  integers, 
i  ]i 

r,  (i  =  l92,.00m)  and  c    (j  =  1  .  ? ..  ,  „,r»)  be  positive 

3-  J 

integers  such  that  £  r .  as  2-  c 

Determine  values  of   x  ,   which  minimize  the  sum  /_  x  A,  ,      subject  to 

x   >  0 

ij  ' 

J?-  X  ,  „  a  C  . 
<     lj      J 

2_    x.  =  r . 

7"   U 


60 


s  pape 

We  work  with  the  cost  -  "•   x  D  -  L  .  ,  J 

problem;,  we  will  disc  hie  ma  tit I     tiling 

red  lin     ind  we  wi  i  1 
of  the  matrix  by  means  of  astei  In  add       we 

will        to  eacl 

which  may  be  changed  t  J 

of  the  matrix  who  s sent la 1  ( 

elements  will  always  be  zei  .he 

number  c  -  Z-  x,   \  jth 

column  atJthat  stage s  and  t':  'JT    '-    ■      is  *-he  discrepancy 

of  the  ith  row.   These  discre--  <e; 

when  they  all  vanish, 
problem, 

tall  construct  a  diagonal  >  rix  and  place 

the  d  „        and  the  corresponding  x,  ,     a*:,  gonal.   Initial 

discrepancies  (the  r,  and  c„)  will  b«  :t  and  above  the 

matrix.   Discrepancies  at  st  [thin  dotted 

lines  with  row  discrepancies  to  crepancies  below 

the  cost  matrtXc                -.ow  <  A  . ;  column  dis- 
crepant Lesa  as  ^ 

clutter  zero  quotas  will  not  be  show.  uplied  by 

a  blank  above  the  diagonal  in  each  cell  where  a  •^  i  i  is  assigned. 
Preliminaries 

All  quotas  are  zero9  no  lines  are  cove red ^  are  starred  or 
primed o 

Subtract  from  each  row  of  the  matrix  D  its  smallest  entry.,  then 

subtract  from  each  column  of  the  resulting      t  st  entry. 

Find  a  zero  Z  in  the.  matrix.   If  the  discre  of  both  its 

row  and  its  column  are  positive,  increase  the  quota  assigned  to  Z 

smaller  of  the  two  discrepas  ■  Rgpeat 
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for  each  zero  In  the  matrix,,   Then  cover  ever-  discrepancy 

Step  lo  Chocs*      --covered  zero  and  prime  ■•  ,  le  row      Ln- 

ing  it.   if  the  discrepancy  of  this  row  is  .  2. 

Otherwise  (if  the  tow  discrepancies  Is  zex~  the  row.   Then  scar  each 

-  —  ntial  zero  Z  in  I  Vs   column 0   Re- 

peat seres  are  covered.  Go  ti 

Step  2.   There  is  a  se  :arred  and  primed  zeros  con- 

icted  as  i        Let   Z   denote  t  Ly  one), 

o 

the  0*  in  Z  Bs  colun  Let   Z„  denote  til 

o 

Z   denote  the  0*  in  Z?°s  s      Similarly 

com  ii     atil     sequence  stops  at  a  0°,  Z~  ,    I  -h  has  n  its 

column.   Since  no  column  contains  more  than  oi  and  no  row  more  than 

one  0°s  this  sequence  is  unique .   It  ms :.  Lement.   The 

discrepancies  of  Z  "s      s  posit ivt  Zy,  '  's  column  is  positi 

and  t      •     signed  each  0*  of  the  Z              Z~.    is  positive. 

h  be         Lest  of  these  positive  nu  -  Increase  the  quota  of 

each  0"  In  the  seqv     by  h8  and  decrease  0*  In  t 

rase  all  asterisks  and  primes.   Un<  sl'i  rows  and 

cover  every  column  whose  discrepancy  is  zex  1. 

Step  3.   Let  k  den        smallest  n  :he  matr; 
will  be  positive.  Add  k  to  every  covered  row  and  subtract  k  from  every 

uncovered  column.   Retur     "tep  ls  u  -my  asterisks,  prin 
or  covered  lines. 

Our  sample  problem  will  be  stated  ij  i  and  solved 
according  :            hm. 
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:  30 

4 
c   =  40 
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Results  of  "preliminar 
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4/6 


5" 

5" 

3 

Z 

5 

4     '+ 

.' 

l 

Q 

z 

3  i  -ne  s , 


k  h  r 


i. 


20 


/ 


-/     '  7o 


/o 

,;  o\    i       z 

.:.A^o.  g 

> 


O  ■> 


4- 

/ 
? 


40 


Ac 


o 


4J 


0 


3a 

/  o    -  /   I 


o\  o\&o\  c  :  o  ; 

— i — i — i — i — i 
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of  Step   1 


V, 


10 


40 


10 

o 


70     30 


Ac 


4  .    q        0 


■  3 

2o 

.   o     i  !    z     i      o,  o 


10/  : 

,3  d-o.y-.  ^i 


^  :  o 


i  i 


C     (So  :     O  :   O    ! 


Step  t  _,  ?.         (3,3) 

30 


r\    Zo  ■  4-0 


/a 


10 


7o 


- 


U- 


i-0 


JO 


<kc 


.,_ .  — , 


,2a, '  fa.  - 


,'1«4<',|  /I    -' 

/3w  o  /O'/ZyZ^ 


\,3\/oVzyoVt 


&t 


0 


0    ' 

.-:       * r 


_0 
30 


I  t 


I  I 


*j    :  O  :  0\3o:  o\  o: 


Our  i mg   the 


■ 


tep  3  oui        20   Applying      3  chan>: 

ent 

step 

2  rfith  h  : ion  is i 


10 

Zo 

4-0 

\0  \ 

i 

30 

3o 

30 

. 

The  cost  is  430 s  as  b«  • 

The  simi.larir.ies  between  the  latter  two.  algorithms  is  apparent  in 
the  manner  in  whic'         : ions  developed 0 


Ml 
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Hill  I         Enco,  New  York  ,  19! 
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APPEN; 
SOLUTIONS  OF  A  MAXIMAL  FLOW  PB 

lc   General. 

In  this  appendix  we  shall  present  --  small  caps      ted  netwoi 
and     rmine  the  maximal  flow  therein  by  using  the  techniques  described 
in  Chapter  III.   We        -  rk  the  problem  first,  by  ehe  flooding  tec 

and  second  by  using  the  minimal-cut  procedure,,   In  each  instance,,  we  shall 
restate  the  essential  features  of  the  te      ae  employe 
20   The  Problem. 


For  the  network  N  shown  above  determine  the  maximum  st?ady- state  flow, 

from  origin  a  to  terminal  b0 

3  ding  Techn 

•rigins  as  Lc  flows  to  all  arcs 

leaving  the  origin  to  saturate  tl  Lon  points  as 

new  origins  am!  greatest  capacity 
constraint,  schedu         whenever  possible  in  the  following  order. 

1„   "Forwar              :ticn  outgoing  arcs, 

°.    "Lateral'        -.it her  points  st. 


3,    'BottleneckecT'  -  If  units  are  left  over  at         1  and 
E  all  outgoing  and  lateral  .rated. 

Continue  the  above,  procedure 
complete  network  and  reaches  the  tern 

Eliminate  bottlenecks  by  returning  a  I         nil   to  the  origin. 
The  validity  of  the  solution  can  be  che<      /  insp>        If  a  maximal 
flow  has  been  achieved,  there  will  be  nc  unsaturated  path  ex- 

tending from  the  origin  to  the  terminal, 

'  lyi  -  ff ri  s  Flooding  Th 


2'- 


lo 


> 


dL"io 


9 


i 


f      r 


J- 


r\ 


3' 


lie 


'4 


to 


k 


b 


The  above  represents  the  res  first  two  stages  of  flood ir\ 

Bottlenecks;1   6  at  junction  >         jun<  Lon  d. 

All  arcs  arriving  at   points  two  arcs  removed  fi     e  origin  are  saturated, 
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;her  application  yields  the  ted 

arcs  is  shown 


R. 


0 

'5 

\% 

<k 

1 

J 

One  additional  unit  is  bottlenecked  at  point  f„ 

Return  bottlenecked  units  to  origin,,  and  the  maximum  flow  is  seen  to 

be  28  units, 

5,    Tb  it  Procedure o 

Let  1  be  the  chain  joining  a  and  b  which  is  topi  [mpose 

as  lar<><   i  flow  as  possible   (T,  k,)  on  this  cl      thereby  saturating 
one  of  Lts  arcs. 

Subtract  k  from  the  capacity  of  each  arc  of  To 

i 

Remove  the  previously  saturated  arc  whose  capacity  is  now  zero.   Record 

k  .   Continue  this  procedure „   Eventual lys  the  graph  will  disconnect  and 


the  maximal  flow  is  established  as; 
F  b  X  K 


i  =  1, 


where-   k   is  the  amount  of  flow  necessary  to  saturate  the  topmost 

th 
chain  in  the   i   step,  and  n  is  the  number  of  steps  required 

for  the  graph  to  disconnect. 
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6.    Solution  by  the  Minimum- Cut  Procedure 
Step  1 


T  =  (a,c8e,h,bs 
Step  2 


a 


n 


Zd 


EX 


k1  -  8 


T  *   (asc,f,e,hsb . 


k2  =  l 
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o   3 


JL 


a 


(asc,f„j 
Arcs  ej  and  hb  can  no  1c      antribmte  to      hain  flow  from  a  to  b, 

Step 


C 


It 


a. 


5" 


c!L 


* 


7 


» 


-  <a,c,d,f,j 


Step  6. 


C 


2 


CL 


-Xo~~ 


*— 


1  9 


/0 


T     ^    (a,c,ds  E 
Step   7 


,2 


CL 


fc 


^6 


T.    =    (a,d,f,g,j,b) 
Step  8 


Tg   -    (asdsg9j,b) 


- » 


/o 
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Step  9= 


a. 


Ho 


I 


-I 

3 

/ 

<o 

/ 

■'  10 

? 


L 


ft 


Tg  =  (a,drgsk8j. 

Step  10o 


A. 


Z 


5 


y 


10 


T1Q  =  (a,dsg,k,b) 


\         ~    1 


On  completion  of  thi         s  g^aph  disconmei 

The  capacity  of  N   is  given  as:' 

Cap(N)  =  5k,   =  8+1+3+2+3+2+2+2+3+2  -  28  units. 
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APPENDIX  C 

GRAPH  THEORY 

In  Chapters  II,  III,,  and  IV  of  this  thesis  we  have  used  the  language 

of  Che  mathematical  theory  of  graphs.   We  include  here  defi     >ns  and 

theorems  from  this  theory s  in  order  to  amplify  and  r.larify  some  of  : 

statements  in  the  text. 

1.   DEFINITIONS. 

A  graph        a  finite^  one-dimensional  complex,  consisting  of 
vertices  af,bj,c, . . .  »e  and  arcs  c^(ab),  #»ac)8„..8  cf(ce)0 

An  arc  c^  (ab)   joins  it  end  vertices  a  and  b;  it  passes  tht 
no  other  vertices  and  intersects  no  other  arcs  except  in  verticc 

A  chair,  (path)  is  a  set  of  distinct  arcs  of  G  which  can  be  arranged 
as  c^(ab)s  /S  (bc)}...,  cT(gh)  where  the  vertices  asb8c8,„,  distinct. 
A  chain  does  not  intersect  itself  and  it  joins  its  end  vertices,,  a  and  h. 

A  cycle  is  a  set  of  distinct  arcs  that  can  be  ordered  as  cC(ab)8 
P   (bc),...„  y(ef)^  cf (£&) »    the  vertices  being  distinct  as  in  the  case 
of  a  chain, 

A  graph  is  connected  if  each  pair  of  vertices  is  joined  by 

A  forest  is  a  graph  containing  no  cycless  and  a  tree  is  a  connected 
forest. 

We  may  associate  with  each  arc  of  a  graph  a  positive  number  called 

;raph  G  together  with  the  capacities  of  its  individual  ar< 
■ twork, 

ted  or  directed  graph  is  a  graph  in  which  each  arc  has  a 

ed.   In  such  a  graph,  C^   (ab)  and  (X  (ba | ,    they  both 
ex i  tinct. 

some  applications  with  oriented  graphs  it  is  desirable  to  dis- 
tingoi      rain  subsets  of  vertices  as  origins  and  certain  others  as 
desl  names  used  in  this  connection  are  or-         ter- 

rid  sources  and  s^nks.   The  word  terminals  is  also  often  us 

ilnguish  a  subset  ©f  nodes  in  a  network,,  the  remaining  nodes  being 
;.  ion  points. 
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A  two-terminal  network  is  a  network  in  k  r.  ?d 

and  distinguished  both  from  each  othe i 

e  point 8  is  i 
called  the.  sink. 

A  chain  flow  from  a  source  a  to  sink  b  is  a  c 
of  a  chain  joining  a  and  b  and  a  non-nej 
flow  along  C  from  a  to  b„ 

A  network- flow  is  a  collection  of  chain  flows  which  has  the  pr> 
that  the  sum  of  the.  numbers  k  of  all  chain  flows  that  contain  any  arc 
no  greater  than  the  capacity  of  that  arc.   If  equality 
said  to  be  saturated  by  the  flow, 

le  flow  in  a  two-terminal  network  is  said  to  be  uni- directional  i 
among  all  the  chains  from  source  a  to  sink  b  ail  chai 
same  direction  in  any  one  arc. 

1    DISCUSSION 

It  is  often  desirable  to  associate  &n   abstract  g 

logical  graph  QJ  ,   This  can  be  done  b\ 

a  point  in  Euclidean  three  space  and   each  arc  as  a  .  ts 

end  vertices:   In  keeping  with  the  definitions  ab  < 

intersect  except  at  vertices.   A  graph  1     Ld  to  be  pi 

1  -  i 

ed  topological  graph  is  planar.   Whitney  UJ  gives  other  criteri 

determining  whether  a  graph  is  planar. 

following  two  theorems. 

Theorem  1.   A  graph  is  planar  ifs  and  ow.lv  if,  it  has  a  dual. 
Theorem  2.   A  graph  is  planar  if,  and  only  if,  it  contains  neither  : 
following  two  graphs  as  subgraphs: 

(A)  G'1      consisting  of  five  nodes  a,b,c,d8e,  and  arcs  connecting 
of  these  lodes  to  every  other  one  by  a  single  arc  or  a  s- 

. chain.  (Note.  A  suspended  chain  is  a  chain  in  wh:  ad 

vertices  are  met  by  more  than  two  arc 

(B)  G9  e  consisting  of  two  sets  of  three  nodes  each  a,l     as  : 

and  arcs  connecting  each  of  the  nodes  of  the  first  set  to  everj 

of  the  second  set  by  a  single  arc  or  a  .     nded  chain. 

1,   A  topological  graph  is  planar         \  be  snappy  - . 
a  sphere. 


Graph       a  rep- 
graph  theory  - 

Representations,, 

An  obvious  repress      >n  of  a  graph 

are  used  to  represent  nodes  and  lines  tc 

be  shown  as  numbers  adjacent  to  the  arci 

tion  is  usually  denoted  by  arrows,, 

Another  representation  of  a  gi 
the  rows  and  the  columns  represent' 

graph  G  by  a  matrix  L  =  Ji  Jt  S 

in  the  graph  and  is  zero  otherwise    I     Ld  r 
capacities  by  replacing  the  lus  with  number 
In  an  oriented  graph    Lj  may  differ  ftorc  JL  j 
positive  and  the  other  ze  i 

Another  matrix  represent.,'     of  a  grap« 
this  matrix  the  nodes  are  represented  by 
say  an  arc  is  incident  on     de  (and  vice  - 
the  node,   In  an  oriented  (directed;*  graph  we  say   that  an  - 
in  to  a  node  is  positively  incident  on  thafl 
from  a  node  is  negatively  incident  on  it„   3 
cidence  by  a  +1,  negative  i   i  -nee  by  a  -1^,  an  ■ 
then  in  the  incidence  matrix  of  a  directed  grap 
+  1  and  one  -1,   Rows  which  represei  • 
row  of  exclusive]/  +1°  s  would  indicat     \rminal» 
be  represented  by  rows  containing  l's  of  mi: 


BIBLXOGRA 

lo   Whitney s  H.*  "Non-Seperable  and  Planar  Graj 

the  American  Mathematical  Society,,  Volume  34s  pp.  .- . 


■  1        1 ... 


IhesC79 

A  review  of  techniques  in  transportation 


III  ill!  II II  III  111  III  fllSlffra 

3  2768  001  02451  6          fiHSilSifiilinMliiiwI^^^^^ffi 

DUDLEY  KNOX  LIBRARY               I'tt'lWl  IM  \v'l\\  M^"^''^^  *  J  *  I  *1 1  i  *  *  -■  ill*l  V  \^*  -  *  t^  i^lVlf  ill« 

I 


1  litliuiMiffliiii''---  ■ 

ftl.l IP     1   11  ill  tiiiii 

i  I  I     1     a    Hi 


!   '  il    II  '    fJli! 


H  till   J  i  ill  1  ilill'l  11  IIIMhr  ItrIiIILI' 
rlHl  llPiPil!  11  1  111     nffl'll mfi    I 

'I  !!'•  it    il       11  1  ill 

Pi    M  li  '  1  Mil  »  1  It  1 1  In  111  1    hIIIHhuihI 

i     I         1          11 

I'M  I1 1.  il>{  j  lii!  1 1  Ml  ]l  (t'Si'i  i']Hii}iliu»n  IliiH 

1             '    '     ' ■             '  ' '  • ■     ' ■  1  '■ '  '  '  ' '  l'i             !  '  ' 
i'i'  m'i   llUnllU'lli Ml!  1  iill'l'I'M'I'l  i'JiIh'iIiiiiiI'I' 

H^HH  i'l'lil'lili  i  1  i     ('  ■  I'  1j  n!  ni'l 

1  i  i  i< ! . :  1 1 \\  i  i '  i  llfi'i  lis  i    iUilHltiii     1 

